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Abstract 

We study the 0{N) linear sigma model in 1+1 dimensions by using the 2PI formalism 
of Cornwall, Jackiw and Tomboulis in order to evaluate the effective potential at 
finite temperature. At next-to- leading order in a 1/N expansion one has to include 
the sums over "necklace"' and generalized "sunset" diagrams. We find that — in 
contrast to the Hartree approximation — there is no spontaneous symmetry breaking 
in this approximation, as to be expected for the exact theory. The effective potential 
becomes convex throughout for all parameter sets which include N = 4, 10, 100, 
couplings A = 0.1,0.5 and temperatures between 0.3 and 1 (in arbitrary units). The 
Green's functions obtained by solving the Schwinger-Dyson equations are enhanced 
in the infrared region. We also compare the effective potential as a function of the 
external field (f) with those obtained in the 1PI and 2PPI expansions. 
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1 Introduction 



The O(N) linear sigma model at finite temperature has a long-standing history, in par- 
ticular as a basic model for a quantum field theory with spontaneous symmetry breaking 
[H El El E] Early investigations beyond the classical level have been based on includ- 
ing one-loop quantum and thermal corrections. These studies have been centered around 
the discussion of the one-loop effective potential Kff(0) where is the mean value of the 
quantum field in a sense being defined more precisely by the effective action formalism, 
summing up one-particle irreducible (1PI) graphs. A next class of approximations include 
bubble resummations, usually obtained from the 1PI formalism by introducing an auxil- 
iary field. Within such a framework the effective potential has been computed to leading 
j2 El E] and next-to- leading order (NLO) of a 1/N expansion [BJ. At leading order in 
1/N one does not find spontaneous symmetry breaking in 1 + 1 dimensions, while in 3 + 1 
dimensions the effective potential is flat for \<f>\ < v |I| 3 . At next-to leading order there 
is no spontaneous symmetry breaking in 1 + 1 dimensions either. These results correctly 
reproduce a general property of theories in 1 + 1 dimensions, where spontaneous symmetry 
breaking is prevented by the non-existence of Goldstone bosons [7j . 

Another approximation including some nonleading terms of a 1/N expansion is the 
Hartree approximation. It is usually motivated by a self-consistency of one-loop quantum 
corrections. It has been studied in 3 + 1 dimensions by various authors, in thermal equilib- 
rium [HJ El Ell HH H2] , and out-of-equilibrium [T3|. The model with spontaneous symmetry 
breaking is found to have a phase transition of first order towards the symmetric phase 
at high temperature. In 1 + 1 dimensions the Hartree approximation displays sponta- 
neous symmetry breaking and symmetry restoration as in 3 + 1 dimensions, an obviously 
unphysical feature. 

If one wants to go beyond the large- N and Hartree approximations there is a variety of 
choices. One of those is the use of the two-particle irreducible (2PI) formalism of Cornwall, 
Jackiw and Tomboulis (CJT) ^3J within which one may include higher loop corrections, 
higher orders in 1/N etc., as in a 1PI expansion. But furthermore it becomes necessary to 
solve Schwinger-Dyson equations for the propagators in order to sum up certain classes of 
Feynman graphs. Apart from the fact that this may be technically demanding, there still 
is the obstacle that for 3 + 1 dimensions renormalization has been discussed up to now only 
for vanishing external fields [T^l EH1 EH EH] ; this obstructs the discussion of spontaneous 
symmetry breaking, which here is our main interest. As here we discuss the model in 
1 + 1 dimensions, the problem of renormalization does not arise and explicit numerical 
computations can be performed. This is indeed the main subject of this publication. 

A technically less demanding approach is the so-called two-particle point-irreducible 
(2PPI) resummation introduced by Coppens and Verschelde ^Hl E01- Here, instead of 
treating the Green's functions as variational parameters, one just introduces variational 
masses like in the Hartree approximation. This implies that the resummation is only over 
local insertions, the 2-particle point-reducible graphs, i.e., graphs that fall apart if two 

3 This follows from an end point extremum of r(</>, M 2 ) at M 2 = 0. 
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lines meeting at the same point (the 2PPR point) are cut. This approach is identical to 
the Hartree approximation if only one-loop 2PPI graphs are included. For this formalism 
renormalization has been fully discussed |2l] : symmetry restoration at finite temperature 
has been investigated in 3+1 dimensions in an approximation including the sunset diagram, 
for the case N — 1 j22] and for the O(N) model with arbitrary N |23j . 

The 2PI formalism has the advantage that at a given order of the loop expansion it 
resums a larger class of Feynman diagrams than the 2PPI or Root's auxiliary field scheme. 
From a variational point of view it allows for a more "flexible" propagator whose self-energy 
insertions are nonlocal in space-time and thereby become momentum dependent. 

It is the purpose of the present work to study the 2PI formalism at next-to-leading 
order of the 1/N expansion, an approximation usually denoted as 2PI-NLO for short. For 
the O(N) model the next-to-leading order 1PI diagrams have been identified in the work 
of Root [H] using an auxiliary field method. They are depicted in Figs. ^ and |H1 in a way 
that graphically differs from Root's presentation. The 2PI-NLO diagrams are obtained by 
omitting those which are two particle reducible, they have been discussed in Refs. (2H EE] ■ 
The sums over all "necklace" (cf. FigJTJ and "generalized sunset" (cf. Fig. |2J) diagrams 
can be done in closed form. 

We have evaluated the effective potential at finite temperature in the 2PI-NLO scheme 
by solving the Schwinger- Dyson gap equations. For comparison we will also show results 
obtained in the other schemes introduced above. 

The consideration of <3> 4 theory in 1 + 1 dimensions has recently received some interest 
in the context of nonequilibrium dynamics [2H I2E1 123 EE I2H]- A comparison [23 EH] 
of various approximations for the 0(1) model, i.e., a model with a simple double-well 
potential, has displayed some still not well- understood features. The 0(1) model has 
thermal kink transitions [HU 1221 EH1 EH] that lead to a symmetric phase in the exact 
theory. These are out of reach within a perturbative treatment in thermal equilibrium, 
whereas they may play (or seem to play) a non-perturbative role in nonequilibrium physics 
even in the approximations considered here. Our present calculations are not directly 
relevant to this case. However, the comparison of the various schemes for N > 1 should 
be of interest as well, both in equilibrium and out of equilibrium. Here we provide the 
equilibrium computations, in expectation of their nonequilibrium counterparts. 

The plan of the paper is as follows: in Section E] we present the general formulation of 
the model and of the 2PI formalism at NLO-l/iV. We also present all the formulae used in 
computing the finite-temperature effective potential in the NLO-l/iV approximation. This 
discussion includes the Hartree approximation. The analogous derivations for the 2PPI 
and 1PI (Root) scheme are presented in Appendix 1X1 and IB1 In Section E] we discuss our 
numerical results. We end with some conclusions and an outlook in Section 0] 
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2 Basic equations 



The Lagrange density of the linear sigma is given by 

C=±dp$.P*-±(&-Nv 2 ) 2 . (2.1) 

The 2PI effective action formalism introduces two variational functions: an external field 
4>(x) and a Green's function Q(x,x'), which arise as Legendre conjugate variables related 
to a local source term J(x)Q(x) and a bilocal source term K(x,x')<&(x)<&(x'). Assuming 
that the expectation value of $ 2 will scale as N we introduce the external field <fi with a 
scale factor y/N, i.e., (<&) = y/N(f>. We then obtain the classical potential 

W0) = ^ (N<P 2 - Nv 2 ) 2 = N± (0 2 - v 2 ) 2 . (2.2) 

The one-loop part of the effective action is given by 

IVioop[0,£] = '-Tilng-'/g, 1 + ^Tr [iV- x Q - 1] . (2.3) 

Here Q, Qq and V are matrices. If we separate the fields in a basis parallel and orthogonal 
to the direction n of the classical field ($) c i = y/Nncj) these matrices become diagonal, with 
one entry for the parallel component and (N — 1) identical entries for the orthogonal ones. 
We denote them as a and ir in reference to the linear sigma model, which presents one of 
the possible applications. With this notation we have 

iV' 1 = -d»d» - A(/,0 2 - v 2 ) . (2.4) 

This is the Klein-Gordon operator in the external field; fj = 3 for j = a and fj = 1 for 
j = ir. We have normalized the Green's functions Qj with respect to free Green's functions 

Gil = -d^ - m 2 fl . (2.5) 

The higher terms in the effective action are all two-particle irreducible vacuum Feynman 
graphs with external lines y/N(J)(x) and with Green's function Q as internal lines. 

In the following we will consider the system at finite temperature. So we perform the 
Wick rotation and replace the integration over the Euclidean p° by the summation over 
the Matsubara frequencies p° = u n = 2mrk B T. Furthermore we consider the effective 
potential instead of the effective action. Then the one-loop term for one species, j = a or 
j = Ti takes the form 

„ 1 T V- r dq f, g 2 + a& + m%Q + S ) -(u )m ,g) 

Vi -"°" = 2 T ^L^r — — (2 ' 6) 

- [m 2 fl + Ej-(w m , q) - X(fj(j) 2 - v 2 )] Gj(uj m , q)} . 
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Here 

g^(u n ,p)=p l + ul + ml (2.7) 
is the reference Green's function, and Ylj(u n ,p) is defined as 

Ej(w n ,p) = Gj x {w n ,p) - g~^(uj n ,p) . (2.8) 

The choice of only introduces an additive, though temperature dependent, constant to 
the effective potential. An obvious choice would be to use the bare physical masses, but this 
would imply m 2 = 0, leading immediately to infrared singularities. As we anticipate that 
the symmetry may not be broken in the NLO-l/iV approximation, we prefer a symmetric 
choice: m 2 = Xv 2 for both j = a and j = it. 

The one-loop potential is divergent. The asymptotic behavior of the integrand is 

So the integral of the one-loop effective action can be regulated by subtracting this term 
and by adding it in regularized form. We have 



1 " l00P ' j " 2 T Aj-oo^l ln S 2 + ^ + - 2 ,o 



m=— oo 



■ [m 2 + Ej(w m , g) - A(/j0 2 - v 2 )] Gj(uj m , q) 



- [m% - A(/,0 2 - v 2 )] B j0 (2.10) 
where Bj are the bubble diagrams which we regularize dimensionally as 

(i-e) 



f dq ( ' 

Bjfl = T I (27r) (i- e) ^o(^ m ,g) 

-m — — ' 



m=— oo 

l 



47T 



2 ,1.1 m I.ol , f dq 

1e + In 47T - In 

e fi 2 



+ I ^T, / • ( 2 - n ) 



We define their finite part Z3j,o,fin using the MS prescription, leaving out the terms 2/e — 
7£ + ln(47r). Then the renormalized bubble diagrams become 

Bj=T jr [<?>„,p) - ^,oK,p)] + B 3fiM . (2.12) 



with 

,2 



-1 , rrij fi i /" c?g 



BjAfin = T" 111 2 + / ^ / • ( 2 - 13 ) 

4tt ^. J *™ exppJqt + mlo-l 



For the renormalization scale we choose //| = m| = Xv 2 . In 1 + 1 dimensions the theory 
is renormalized if the Hamiltonian is normal-ordered. This is not a unique prescription 
as the normal-ordering may be done with respect to different masses of the bare quanta 
|H5| IHfij. A shift in these masses introduces a redefinition of the vacuum expectation value 
v , or of the mass parameter in the symmetric theory, as the change in normal ordering of 
the 4 term introduces a term quadratic in the fields. In the present context we have to 
take the renormalization scale identical for the a and 7r modes as otherwise we break the 
symmetry explicitly. Our convention corresponds to normal-ordering with respect to bare 
quanta of mass Xv 2 . 

With these preliminaries the total one-loop effective potential is given by 

^i-ioop = ^i-ioo P , CT + {N- l)Fi_i ooPj7r . (2.14) 

The first nontrivial term is the double-bubble diagram. As we consider here leading 
order and next-to-leading order contributions we have to calculate this diagram with the 
exact combinatorial factors for O(N). It takes the form 

V dh = -^{[B (T + (N-1)B 7T } 2 + 2[B (T + (N-1)B T }} 

= ^ [(N 2 - l)Bl + 2(N - l)B a B« + 3Bl] . (2.15) 

If only these two contributions are considered, we obtain the Hartree approximation, or 
daisy and super-daisy resummation. In writing down this contribution we have replaced 
the divergent bubbles subdiagrams by the finite ones. This requires mass and vacuum 
energy counterterms which here are fixed by the minimal subtraction prescription which 
eventually may be replaced by a precise renormalization condition. 
From 

^Hartree = ^1 —loop 

+ V dh (2.16) 

we obtain the gap equations 

SfK,p) = -ml + X(3(f> 2 -v 2 ) + ±[3B ff + (N-l)B v ] (2.17a) 

S?K,p) = -m 2 n0 + X(<f> 2 - v 2 ) + A [B ff + (JV + 1)5,] (2.17b) 

When taking the functional derivative with respect to Q n we will, here and below, consider 
each pion "individually", disregarding factors (N — 1) which would cancel in the gap 
equation. In this way the second equation is obtained directly, and the sigma and "each" 
pion field are treated on the same footing. 

The large-iV limit is obtained by omitting the sigma contribution entirely. Then the 
first of the gap equations is omitted as well, and the second one only contains the pion 
bubble with NJ (N — 1) replaced by unity. 

The condition for a nontrivial extremum with respect to <fi becomes 

1 8V H 

= =^ + B-» ! ) = (2.18) 
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with the compound bubble diagram 



B = = 1 [3B a + (N — 1)^] . (2.19) 

Again in the large- N limit the sigma contribution is neglected and the prefactor N/ (N — 1) 
of the pion bubble is replaced by unity. To next-to- leading order in a 1/N expansion one 
has, in the absence of external fields, to sum over all diagrams which are traces of powers of 
the fish diagram, the so-called necklace diagrams, shown in Fig. ^ In order to go beyond 
the leading order we define several diagrams: We denote the fish diagrams by fj(w n ,p); 
explicitly 



Tj (u) n ,p) = T ^ J ^ Qj (w m , q)g j (uj n - uj m , p-q) 



(2.20) 



and denote by !F(uj n ,p) the "average" fish 



H"n,p) = ^ [{N- l)^(cu n ,p) +F a (u n ,p)} . (2.21) 
The sum over the necklace without external fields, see Fig. Q is then given by 

m=—oo 

The subtraction of XJ-(u) m ,p) accounts for the fact that we have already included the 
double-bubble diagrams. In the symmetric theory where the a and tt contributions are 
equal, the factor N cancels, so this expression is of order N° which is NLO. 

The same order in 1/N is achieved by replacing one of the a lines of the necklace by 
N(f> 2 . We will call this class of graphs "generalized sunsets", see Fig. |21 for a graphical 
depiction. We introduce the derivative 

sg.i^P) = n t 1^ ^TTAfM) g ^ n - Um > p - q) ■ (2 - 23) 

J m=— oo v ' 

Note that the factor X/N arises from the derivative of the expression XJ-{u n , p) with respect 
to Qj. In terms of 5Af/5Q a we obtain for the generalized sunset diagram, or sum over 
necklace diagrams with two external lines 



The functional derivative contains a factor 1/N so that this expression is of order iV . 
With these definitions the 2PI effective potential up to NLO becomes 

V cS = Klass + ^1-loop + V db + S + M (2.25) 
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Figure 1: Resummation of necklace graphs. 




Figure 2: Generalized sunset graphs. 



In order to write down the gap equations we need to introduce a further functional 
derivative 



—(u n ,p) = N(f) — — 



SAf 



(0,0) 



Explicitly it is given by 

, N A 2 ~ f 



dq 



-1 



2tt [1 + A^(cu m , q)\ 



2 Qa ( W ra , q) Gj ( UJ n ~ ^ra , P ~ Q) 



2i -A.F(u;„,p) 



With these definitions the gap equations become 



A 



= -m 2 + A(30 2 - v 2 ) + - [(N - 1)3* + 3B a 

SAT. , 5S . 
+ — {u n ,p) + —(uj n ,p) 
oy a oy a 

E n (u n ,p) = -ml + \^ 2 -v 2 ) + ^[(N + l)B w + B a ] 

SAf , SS 
oy n dy n 



(2.26) 



(2.27) 



(2.28a) 



(2.28b) 



All expressions in these equations have been given explicitly above. Finally for the partial 
derivative of V e g with respect to we find 



(2.29) 



We have performed the computations for the 2PPI scheme of Verschelde and Coppens as 
well, including necklace and generalized sunset diagrams. The formulae and their derivation 
are rather similar to those of the 2PI scheme. They are given in Appendix El 

3 Numerics and Results 

The equations of the previous section are easily incorporated into a computer code; we 
solve the gap equations by iteration. The range of momenta and Matsubara frequencies 
was restricted to values smaller than 15 — 20, far above the relevant mass scales. All 
numerical integrals and summations are ultraviolet finite by the regularization presented 
in the previous section. 

For the first step at fixed N,v,X,T and <fi we start with Qj = Qj . The iteration is 
monitored by the values of the average bubble diagram 

B=[(N-1)B 7T + B a ]/N . (3.30) 

If the relative change of this global quantity has become smaller than 10 -8 the iteration is 
considered to have reached the solution. 

For temperatures of the order of v this procedure converges well. Problems arise if 
the temperature is low, typically less then v/5 and for typically less than v/2. In these 
regions the pion propagator becomes large at small (u) n ,p) obviously due to a pion pole 
approaching p 2 = from negative (Minkowskian) p 2 . This situation is close to an infrared 
divergence, and small changes of the pion Green's function result in large changes of the 
various integrals. One can extend the domain of convergence by using an underrelaxation, 
defining the nth step of the iteration of Eqs. ()2.28j) by taking the result of the 

previous step multiplied by a factor of a plus the right hand sides of these equations 
multiplied by (1 — a). We have chosen a typically between .5 and .95. 

In order to be sure of the consistency of our search for extrema, both in the analytic 
formulas and in the numerical computations we have not only computed the results for 
dV c f[[(j),Gj]/d(j) but also the potential Kfr[0, Qj] itself. Here in both cases the Green's 
functions Qj were the solutions Qj of the gap equations; so both V e s and dV e s/d(p are 
evaluated at 5V e g/5Q = and the extrema of the potential Kfr[0, Qj] should coincide 
with the zeros of dV e g[<j>,Qj]/d<l). For the computations in the 2PI scheme we do not find 
any nontrivial minima, so all we check here is that TZ((p) has no zeros. For the Hartree 
approximation and for the 2PPI scheme, the extrema of the potential and the zeros of 
7Z(4>) agree within the avalailable accuracy. 

In Figs. H3 and 0] we show the temperature dependence of the 2PI effective potential, 
for N = 4 and A = 0.5 and A = 0.1. The effective potential shows no sign of an inflection 
point. We also plot the Hartree result for one of the temperatures, in the neighborhood of 
the critical temperature (of the Hartree approximation). 

In Fig. 03 we show the comparison between the Hartree, the 2PPI-NLO and Root's 
NLO approximation (denoted as 1PI-NLO for short, cf. Appendix El) and the convergence 
towards the leading order large- N results. These results are for v = 1, X = 0.5 and T = 0.5. 
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Figure 4: 2PI effective potential at NLO-l/iV for iV = 4 and A = 0.1. 
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As one sees, the effective potential in the 2PPI approximation is always in between those 
of the Hartree and 2PI-NLO approximation. The 1PI-NLO effective potential is close to 
the 2PI results. The latter two approximations differ only by terms of next-to-next-to- 
leading order (NNLO); it is surprising, nevertheless, that these terms are small already at 
moderate values of N. They are important for the evolution out of equilibrium, as they 
include scattering of quantum fluctuations. As for finite N one expects differences to the 
large- N result, the deviations of the various approximations from this limit do not a priori 
establish any "ranking"; however, the spontaneous symmetry breaking displayed by the 
Hartree and 2PPI-NLO approximations is certainly unphysical. 
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(c) N = 100 

Figure 5: Comparison of the effective potentials obtaind in the Hartree, 2PPI-NLO, 1PI- 
NLO, 2PI-NLO and leading order large-iV approximations for A = 0.5, T = 0.5 and 
different values of N. 



In Fig.lHlwe display the spectral behavior of the effective momentum-dependent masses 
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ml + T,(u n ,p), see Eqs. (|2.28|) . at the equilibrium point <fi = 0, for various temperatures, 
N = 4 and A = 0.5. We plot these quantities versus the "Euclidean momentum" p E = 
^Ju)n+ p 2 - The curves for various Matsubara frequencies u n are close together, implying 
an approximate rotation symmetry in the Euclidean (u,p) plane. Whereas the large- 
momentum behavior is determined by an effective constant mass given by the tree level 
and bubble contributions, the effective masses at low momenta are considerably smaller, 
so that the various loop integrals are infrared-enhanced. The infrared enhancement can 
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Figure 6: Momentum-dependent masses, cf. Eqs. (J2.28)) . for different temperatures. N = A, 
A = 0.5. Here pe = >/uj% + p 2 . 



also be seen in Fig. where we plot, for A = 0.1, the values of the effective masses at the 
equilibrium point = 0. We display m\ + E(0, 0) and m\ + E(0,p max ), where p max is our 
pragmatic momentum cutoff. The values near p = are much smaller than those at p max , 
which can be considered as the asymptotic ones, see Fig. |H1 



4 Conclusions 

We have analyzed here the O(N) linear sigma model in the 2PI formalism at next-to- 
leading order in a 1/N expansion. In order to do so, we have solved the coupled system of 
Euclidean Schwinger-Dyson equations for the a and 7r propagators. 

We have found that within our range of parameters the 2PI approximation does not 
display spontaneous symmetry breaking, as to be expected in a 1 + 1 dimensional model. 
Furthermore, we find that the thermal propagators are enhanced in the infrared region. 
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Due to momentum-dependent self-energy insertion the effective masses at low Euclidean 
momenta are smaller than those at large momenta. 

We have compared the effective potential as a function of the external field to the 
ones obtained in several other approximations. The 1PI-NLO approximation is close to 
the 2PI-NLO one even for small N, whereas the potential obtained in the 2PPI-NLO 
approximation lies between those obtained in the Hartree and in the 2PI approximations. 
All potentials are above the large- N effective potential and, of course, converge to it at large 
N. The Hartree and 2PPI-NLO approximations display spontaneous symmetry breaking, 
unphysical in 1 + 1 dimensions. This is clearly an artefact of these approximations. 

It would be of interest to extend the analysis of the O(N) linear sigma model in the 2PI- 
NLO approximation to nonequilibrium systems. For the O(N) model out of equilibrium 
the 2PPI approximation (to two loops in the effective action) displays a symmetric phase in 
the late time behavior [SZ|, in contrast to what we have found here in thermal equilibrium. 

Coming back to the simulations in 1 + 1 dimensions for the case N = 1 [2E1 EH] spon- 
taneous symmetry breaking was found in the late time behavior of the Hartree approxi- 
mation, whereas in both the 2PI and the 2PPI the mean field approaches the symmetric 
phase at late times. In contrast to that, the 2PPI approximation in equilibrium displays 
spontaneous symmetry breaking, and so does the 2PI approximation both with correct 
N — 1 combinatorics and with 1/N combinatorics [HB]- This is due to the fact that these 
calculations do not encompass kink transitions which are of order 1/A. So there remain 
some issues to be clarified. 
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In this work we were concerned mainly with the phase structure of the model in the 
2PI scheme at next-to-leading order of the 1/N expansion. It would of course be useful 
to extend this analysis to 3 + 1 dimensions using recently developed renormalizations 
techniques for 2PI-resummed perturbation theory HOI El CHI EEj ■ Furthermore it would 
be interesting to compute various thermodynamic functions and transport coefficients [J01 
I^T| . in view of their relevance for heavy ion collisions. 
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A The 2PPI formalism to NLO of a 1/N expansion 

In the 2PPI formalism ^3120] one introduces a quadratic source term K(x)Q 2 (x) which, in 
contrast to the 2PI formalism, is local. The Legendre transformed variables are then <p(x) 
and A(x). For a comparison between the 2PI and 2PPI formalism see, e.g., Appendix A of 
Ref. j2Bj- For a translation invariant system and A are constants. The effective action 
then becomes an ordinary function of these variables. For the O(N) system, the classical 
field is a vector $ c i = ncj). As for the 2PI formalism we use a basis where the a field is in 
the direction n and the N — 1 pion fields span the orthogonal directions. Then A is a 2 x 2 
matrix which becomes diagonal in this basis, with entries A CT and A n . In contrast to the 
2PI formalism, the Green functions take the simple form 



g j (p,-n)- p2 + u t + M] (A.1) 



with the effective masses 



Ml = A [30 2 - v 2 + 3A a + (N - 1)AJ (A.2a) 
Ml = A [<P 2 - v 2 + A CT + (N + 1) A,] . (A.2b) 



Here the local mass insertions Aj are given by 



np2PPI 

q 



^ = 2 ^hr ■ (A.3) 
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where Pq PPI is the sum of all two particle point irreducible graphs. These are all graphs 
that do not fall apart if two internal lines meeting at one point are cut. They are computed 
with the propagators defined in Eq. (|A.1J) . As for the 2PI formalism we have considered 
the derivatives with respect to Ail as the derivative with respect to the mass of one of the 
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(N — 1) pion species, thus leaving out factors (N — 1) which anyway would cancel in the 
final formulae. 

It is simpler to express the potential in terms of the effective masses M 2 instead of the 
parameters Aj. Then these gap equations can be derived from the effective potential 

UM, Ml, Ml) = NU cl {cf>, Ml, Ml) + U^, Ml, Ml) . (A.4) 

Here U c \ is the "classical" potential [TTH 123] 

Ua = iMir-^'-^^flMlHN-DMl} (A.5) 

[(N + \)M\ + 3(N - l)Mi - 2(N - l)M 2 a Ml + 2N\ 2 v 4 ] , 



lX(N + 2) 



which in this formulation already includes quantum parts through the effective masses. In 
this sense the separation between U c \ and the quantum part U q is superficial. The basic 
contribution to U q is the one-loop contribution 



with 



which explicitly reads 



C/i-ioop = t/f-ioop + (N- l)U?_ loop (A.6) 



^-ioo P = ^Trlog|^. (A.7) 



.7-0 



W -W [dp P 2 + ^ n + Mj 



In regularized and renormalized form it becomes 

p 2 + u 2 n + 



p 1 + + m^ 

2 ™2 



Mj-mj fi \ 1 



p 2 + wj* _|_ m ^ 



} + \{M 2 - mlM . (A.9) 



The regularized bubble integrals have been defined in Eq. (j2.11j) . If we include just t^i-ioop 
we again obtain the Hartree approximation with Aj = Bj. The double-bubble diagrams 
are included here in the U c \ via the effective masses. 

Going beyond the Hartree approximation in the 2PPI formalism in a strict 1/N expan- 
sion we only have to take into account "necklaces" and "generalized sunsets" and a 2PPI 
(but 2PR) combination of them. The necklace Af takes the same form as for the 2PI for- 
malism, see Eq. (I2.22|) and FigHJ where of course the fish diagram, Eq. (j2.21|) . is computed 
with the propagators of Eq. (jA.lj) . The generalized sunset contribution is replaced by a 
more complex set of graphs. As discussed by Root 0, at next-to- leading order of 1/N 
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the graphs with external lines have the form presented in Fig. |HJ with alternating sigma 
propagators and necklaces. As the resemblance to a sunset becomes now very remote we 
refer to them as chain diagrams. These are summed up in the form 



C = \ T E / ^Ml + S( Wm ,g)] 

m=— oo 

where S(u m ,q) is the kernel of the sunset diagram, Eq. (|2.24J) : 



(A.10) 



S(u m ,q) = -2X(f) 



1 + \F(u) m ,qy 



(A.ll) 



In the 2PI formalism the higher powers of insertions S are included automatically, here 
the summation has to be done explicitly. 







Figure 8: 2PPI or 1PI "chain" contributions to the effective action. 

If the graphs defined in the previous paragraph are introduced the quantum part of the 
effective potential takes the form 



U q = t/i-ioop + M + C 
and the insertions Aj in the gap equations are given by 

DM dC 



dM) dM) ' 



(A.12) 



(A.13) 



Explicitly the new contributions are given by 



dM _ \^ T sr- f dp -^(Q^n) 
dM) ~ 2i\T 



2tt 1 + \F(q, u n ) dM) 



and 



dC 
dM) 



(q,^n) 



dTj{oo m , q) 



(A.14) 



M 2 T E / 2 ! 1 + S(u; m , q) { N [1 + \F(co m , q)f dM) " Ga{Un > q) 



\F{u m ,q) 2 

-GA^m,q) 



\J~(uj m , q) 



(A.15) 
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Finally T = {T a + (N - 1)^)/N as before, and 

^(p,w„) = -2TY [ p. j . (A.16) 

dM S m J 271 (<? 2 + < + MJ) [(P - ?) 2 + K - w TO ) 2 + MJ] 

The partial derivative of Kff with respect to (ft, divided by Ncj), in analogy to Eq. ()2.29|) is 
given by 

™ = M *-^ + -k% (A - 17) 



with 

N<f> 



777^: = _2 ^? T Z^ / rr^? ; vF? a.is 



B The 2PPI formalism to NLO of a 1/iV expansion 

In a seminal article [H] Root has discusssed the extension to the next-to-leading order in a 
1/N expansion within the 1PI formalism. He introduces an auxiliary field \ which in the 
language used throughout our work can be identified with a self-consistent pion mass .M 2 . 
The classical potential then takes the form 

f/ cl (0,^ 2 ) = ivQ^ 2 (0 2 -, 2 )-^ , (B.l) 

it can be obtained by taking the large- N limit of the analogous potential of the 2PPI 
formalism, Eq. (|A.5|) . when the mass of the a propagator is given by 

Ml = Ml + 2A0 2 . (B.2) 

To leading order in 1/N this would follow in the 2PPI formalism when taking only the 
pions into account in the one-loop terms. In Root's work this relation is kept fixed for all 
orders of 1/N, however. In a strict 1/N counting one may indeed neglect corrections of 
NLO to the propagator, as actually suggested in Root's work by using the leading order gap 
equation. As in the 2PPI formalism the propagators take the tree level form of Eq. (|A.1J) . 

In the following we will rewrite the various terms included by Root in a way that 
makes the correspondence to the other formalisms more transparent. Including the next- 
to-leading order corrections the one-loop term takes the form of Eqs. ()A.7|) - (jA.9|) . with Ml 
given by Eq. (jB.2|) . As in the 2PPI formalism the next terms are the necklace and the 
chain diagrams. In the necklace part only the pion fish diagram is taken into account. The 
definition of the necklace contribution is then identical to the one in the 2PI formalism, 
Eq. f!2.22|) with the modification that J- ' = and that the propagators in the fish graph, 
Eq. f!2.21|) are given byEq. l[P|) . 

The chain diagram C takes the same form as in the 2PPI formalism, see Eqs. (jA.10)) 
and (jA.llJ) . of course now with T = J-" n . In taking the derivatives the sigma propagator is 
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now considered as a function of and 0. The modifications with respect to the 2PPI 
formalism are obvious. We have 



dC 



-\rh 2 r\^ f dq 1 f A dFjju^g) 

A<P 4- 7 2vr 1 + S(o; m , g) \ iV [1 + A^, q)f dM\ q) 

GZ(u m ,q)\ (B.3) 



1 + A J"(cj m , g) 
and 
1 <9C 



N(f) d(p 

Finally the gap equation is 



4 T E / u?w q) M "™*) [l-2A^ CT (^,g) 

(B.4) 



M ^ X ^^ + 2 m + 2 m) (R5) 
where ON jdJ^A\ is given by Eq. (|A.14j) . and the quantity 1Z is given by 
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